Search for: All records

The majority of computer algebra systems (CAS) support symbolic integration using a combination of heuristic algebraic and rule-based (integration table) methods. In this paper, we present a hybrid (symbolic-numeric) method to calculate the indefinite integrals of univariate expressions. Our method is broadly similar to the Risch-Norman algorithm. The primary motivation for this work is to add symbolic integration functionality to a modern CAS (the symbolic manipulation packages of SciML, the Scientific Machine Learning ecosystem of the Julia programming language), which is designed for numerical and machine learning applications. The symbolic part of our method is based on the combination of candidate terms generation (ansatz generation using a methodology borrowed from the Homotopy operators theory) combined with rule-based expression transformations provided by the underlying CAS. The numeric part uses sparse regression, a component of the Sparse Identification of Nonlinear Dynamics (SINDy) technique, to find the coefficients of the candidate terms. We show that this system can solve a large variety of common integration problems using only a few dozen basic integration rules. 

Abstract BackgroundRepolarization alternans, defined as period-2 oscillation in the repolarization phase of the action potentials, is one of the cornerstones of cardiac electrophysiology as it provides a mechanistic link between cellular dynamics and ventricular fibrillation (VF). Theoretically, higher-order periodicities (e.g., period-4, period-8,…) are expected but have very limited experimental evidence. MethodsWe studied explanted human hearts, obtained from the recipients of heart transplantation at the time of surgery, using optical mapping technique with transmembrane voltage-sensitive fluorescent dyes. The hearts were stimulated at an increasing rate until VF was induced. The signals recorded from the right ventricle endocardial surface just before the induction of VF and in the presence of 1:1 conduction were processed using the Principal Component Analysis and a combinatorial algorithm to detect and quantify higher-order dynamics. ResultsA prominent and statistically significant 1:4 peak (corresponding to period-4 dynamics) was seen in three of the six studied hearts. Local analysis revealed the spatiotemporal distribution of higher-order periods. Period-4 was localized to temporally stable islands. Higher-order oscillations (period-5, 6, and 8) were transient and primarily occurred in arcs parallel to the activation isochrones. DiscussionWe present evidence of higher-order periodicities and the co-existence of such regions with stable non-chaotic areas in ex-vivo human hearts before VF induction. This result is consistent with the period-doubling route to chaos as a possible mechanism of VF initiation, which complements the concordant to discordant alternans mechanism. The presence of higher-order regions may act as niduses of instability that can degenerate into chaotic fibrillation. 

Conduction velocity (CV) slowing is associated with atrial fibrillation (AF) and reentrant ventricular tachycardia (VT). Clinical electroanatomical mapping systems used to localize AF or VT sources as ablation targets remain limited by the number of measuring electrodes and signal processing methods to generate high-density local activation time (LAT) and CV maps of heterogeneous atrial or trabeculated ventricular endocardium. The morphology and amplitude of bipolar electrograms depend on the direction of propagating electrical wavefront, making identification of low-amplitude signal sources commonly associated with fibrotic area difficulty. In comparison, unipolar electrograms are not sensitive to wavefront direction, but measurements are susceptible to distal activity. This study proposes a method for local CV calculation from optical mapping measurements, termed the circle method (CM). The local CV is obtained as a weighted sum of CV values calculated along different chords spanning a circle of predefined radius centered at a CV measurement location. As a distinct maximum in LAT differences is along the chord normal to the propagating wavefront, the method is adaptive to the propagating wavefront direction changes, suitable for electrical conductivity characterization of heterogeneous myocardium. In numerical simulations, CM was validated characterizing modeled ablated areas as zones of distinct CV slowing. Experimentally, CM was used to characterize lesions created by radiofrequency ablation (RFA) on isolated hearts of rats, guinea pig, and explanted human hearts. To infer the depth of RFA-created lesions, excitation light bands of different penetration depths were used, and a beat-to-beat CV difference analysis was performed to identify CV alternans. Despite being limited to laboratory research, studies based on CM with optical mapping may lead to new translational insights into better-guided ablation therapies. 

Abstract AimsThe mechanisms of transition from regular rhythms to ventricular fibrillation (VF) are poorly understood. The concordant to discordant repolarization alternans pathway is extensively studied; however, despite its theoretical centrality, cannot guide ablation. We hypothesize that complex repolarization dynamics, i.e. oscillations in the repolarization phase of action potentials with periods over two of classic alternans, is a marker of electrically unstable substrate, and ablation of these areas has a stabilizing effect and may reduce the risk of VF. To prove the existence of higher-order periodicities in human hearts. Methods and resultsWe performed optical mapping of explanted human hearts obtained from recipients of heart transplantation at the time of surgery. Signals recorded from the right ventricle endocardial surface were processed to detect global and local repolarization dynamics during rapid pacing. A statistically significant global 1:4 peak was seen in three of six hearts. Local (pixel-wise) analysis revealed the spatially heterogeneous distribution of Periods 4, 6, and 8, with the regional presence of periods greater than two in all the hearts. There was no significant correlation between the underlying restitution properties and the period of each pixel. ConclusionWe present evidence of complex higher-order periodicities and the co-existence of such regions with stable non-chaotic areas in ex vivo human hearts. We infer that the oscillation of the calcium cycling machinery is the primary mechanism of higher-order dynamics. These higher-order regions may act as niduses of instability and may provide targets for substrate-based ablation of VF.